Vertex Corrections and the Korringa Ratio in Strongly Correlated Electron Materials 
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We show that the Korringa ratio, associated with nuclear magnetic resonance in metals, is unity 
if vertex corrections to the dynamic spin susceptibility are negligible, the hyperfine coupling is 
momentum independent, and there exists an energy scale below which the density of states is 
constant. In the absence of vertex corrections we also find a Korringa behaviour for Ti, the nuclear 
spin relaxation rate, i.e., 1/Ti oc T, and a temperature independent Knight shift. These results are 
independent of the form and magnitude of the self-energy (so far as is consistent with neglecting 
vertex corrections) and of the dimensionality of the system. 
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Nuclear magnetic resonance (NMR) spectroscopy is 
a powerful experimental probe of the spin dynamics of 
strongly correlated electron materials. An important 
quantity in NMR experiments on metals is the Korringa 
ratio,"'^'^ /C, which is proportional to the ratio of the nu- 
clear relaxation rate 1/Ti to the square of the Knight 
shift, Ks- The Korringa ratio, and the temperature de- 
pendence of 1/Ti and can provide important insights 
into the electronic and magnetic correlations. In this pa- 
per we show that, under certain conditions, the Korringa 
ratio does not deviate from its non-interacting value even 
in strongly correlated electron systems. 

In the diagrammatic formalism of quantum many-body 
theory, the effects of electronic correlations are described 
by the self energy and the vertex corrections. The self 
energy describes the effect that interactions with virtual 
particles have on the propagation of particles through the 
material. Vertex corrections describe the renormalisa- 
tion of coupling constants due to interactions (the name 
arises because coupling constants appear at the vertices 
in Feynman diagrams).- 

An important consideration in the study of vertex cor- 
rections is Migdal's theorem, which states that the vertex 
corrections due to the electron-phonon interaction are of 
order ^JmjM, where m is the electron mass and M is 
the nuclear mass. The Eliashberg theory of supercon- 
ductivity, which improves on BCS theory by replacing the 
BCS effective pairwise interaction with an explicit treat- 
ment of the electron-phonon interaction, invokes Migdal's 
theorem in order to neglect vertex corrections. There 
is no Migdal's theorem for either the electron-electron 
interaction or the electron-magnon interaction.— There- 
fore, understanding the importance of vertex corrections 
in strongly correlated superconductors is of great impor- 
tance and has been widely debated. Hertz el al? have 
shown that the first order vertex corrections due to para- 
magnons is the same order of magnitude as the bare ver- 
tex. More recently, in the context of the cuprates, argu- 
ments have been presented both for— lii^ and agains t^'^°i^^ 
the thesis that vertex corrections are negligible. Some 
theories of strongly correlated superconductors have at- 



tempted to deal with the vertex corrections. For ex- 
ample, for the one band Hubbard model, the FLEX 
approximatioi>i^ consists of a self-consistent summation 
of bubble and ladder diagrams; the latter are the lowest 
order self-constient contributions to the vertex function. 

The relaxation of the nuclei is governed by their cou- 
pling to their environment, which, in a metal, is the 
conduction electrons fii^ Hence, many of the properties 
measured in NMR experiments depend on the transverse 
dynamic magnetic susceptibility of the electron fluid, 
X (-(pt; *^n)- In Matsubara formahsm, this is given byi^ 



X-+(q, it^„) 



/3 



dTe^""^(T,m_(q,T)m+(-q,0)), (1) 



where = l/fc^T is the inverse temperature, r is the 
imaginary time, w„ are the Matsubara frequencies. 



(2a) 



m+(q,r) = ^ E ^l+q,T W^P.i (2b) 



are the =F components of magnetization, and Tt is the 
(imaginary) time ordering operator. 

The Korringa ratic^'^ is the dimensionless quantity 



K = 



1 



(3) 



where 7^? (7e) is the nuclear (electronic) gyromagnctic 
ratio. The Korringa ratio is unity in a non-interacting 
system with a contact, i.e. momentum independent, hy- 
perfine coupling. The hyperfine coupling is momentum 
independent if there is one atom per unit cell, and this is 
an approximation otherwise.^ Further, Korringa showed 
that in such a system l/TiT and Ks are independent 
of temperature^ These three behaviours are often col- 
lectively referred to as Korringa behaviour. However, 
real materials always exhibit some correlations and the 
Korringa ratio may deviate from unity, or 1/TiT or Ks 
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may be temperature dependent i^ii^ii^ In elemental met- 
als the Korringa ratio is typically between 0.6 and 1 (see 
Ref. [H, pp. 156-7). Further, the Korringa ratio is greater 
(less) than unity if the system is near an antiferromag- 
netic (ferromagnetic) instability.^- There are large devia- 
tions from Korringa behaviour in a wide range of strongly 
correlated electron materials including the cuprates,^* or- 
ganic charge transfer saltSjii^i^^ and the heavy fermion 
materialsj^i On the other hand, many materials and 
model Hamiltonians do show Korringa behaviour; includ- 
ing some strongly correlated ones, such as magnetic im- 
purities described by the Anderson model.^^ Therefore it 
is important to determine what conditions may be suffi- 
cient for Korringa behaviour in a correlated system. 

The relaxation rate and the Knight shift can be written 
in terms of the dynamic spin susceptibility. 
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K, 



lim 



2kf 



-0 ilh^ ^ 



\A{ci)\ 



|A(o)|x'-+(o,o) 
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(4b) 



where A{q) is the hyperfine coupling between the nu- 
clear and electron spins, and x'(q, cj) [x"(q, t^)] is the 
real [imaginary] part of the dynamic susceptibility. 

The focus of this paper is on how electronic correlations 
affect the Korringa ratio. Specifically, we investigate how 
vertex corrections modify the Korringa ratio and what 
implications this has for strongly correlated materials. 
We show that as zero temperature is approached the Ko- 
rringa ratio approaches unity if three conditions are sat- 
isfied: (i) vertex corrections are negligible, (ii) the hyper- 
fine coupling is momentum independent, and (iii) there 
exits a energy scale below which the density of states is 
constant. This result holds independent of the strength 
of the electron-electron interactions that enter the self en- 
ergy and of the dimensionality of the system. We also find 
that, under the same conditions, 1/TiT and Kg are inde- 
pendent of temperature. Hence, non-Korringa behaviour 
must result from vertex corrections, a momentum depen- 
dent hyperfine coupling and/or the absence of an energy 
scale below which the density of states is constant. 

Upon substituting ^ into performing the appro- 
priate Wick contractions on the operators, and Fourier 
transforming into frequency space one finds that 



pT 



X_+(q,icj„) = 




(5) 



p + qi 



2„,2 



2/3 



P,Wr, 



xG{p + q,ipm)G{p,ipr„+iujn), (6) 

where r(q, p, zw^j) [shaded area] is the three point 
vertex function and G(p, ipn) [solid lines] is the full in- 



teracting Green's function given, in the spectral repre- 
sentation, by 



r 

G{P,Wn) = / 

J — ( 



2tt ipn ~ El 
where As{p,Ei) is the spectral function, given by 

-2ImS(p,S) 



{E-Sp- RcS(p, S))2 + (Iml](p, E))^ 



(7) 



(8) 



Ep is the dispersion of the non-interacting system, and 
E(p,ip„) is the self energy. Therefore 



X-+(q, iujn) 



fi le 

2(3 ^ 

p,m 



2-K 2-K 



xr(p + q, zpm;p,zpm -^itjJn) 

^ ^(p + q,£;i)A,(p,£:2) 

{ipm - Ei){ip„i + iUn - E2)' 



(9) 



At this stage we neglect vertex corrections, that is we 
set r(q, iu;„; p, ip„) = 1 for all p, q, p„, and ujn- After 
performing the Matsubara summation and the analytical 
continuation icUn —^Lu + irj, one finds that 



x-+(q,^^) = 



dEi dE2 
~2jx'^ 



A,(p + q, El) 



2 ^ 
p 

xA.(p,£;,)"^^^^)-"^(^^), (10) 

nw + El — E2 + ir] 



where np{E) is the Fermi function. 

Using the well known equality l/(x -I- irf) = V{l/x) — 
iTTS{x), where V{y) denotes the principal value, one finds 



lim 



x':+(q,c.) 



2 ^ 
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dE , , 

— ^(p + q,i?) 



x^(P,i^)(-§|).(ll) 



It then follows from Eq. (|4a|) that 
1 kB\A\^ r°° dE 



TiT 



., , , , dnp 

4.^(^H-^ 



(12) 



where p{E) — As{p, E) is the full interacting density 
of states per spin species and we have assumed a con- 
tact hyperfine coupling, A(q) = A for all q. We now 
specialise to the case where their exists an energy scale, 
/csTo below which p{E) is independent of energy. There 
several situations in which Tg may not exist or may be 
too small to be of practical interest, for example, if the 
Fermi energy is at, or very close to, a van Hove singu- 
larity or if there is a (pseudo)gap at the Fermi energy. 
For r ^ To, Tp, where Tp is the Fermi temperature, Eq. 
p2| simplifies to 



1 keW-pHEp 



TiT 



(13) 
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Note that the right hand side, and therefore 1/TiT, is 
independent of temperature. 

For H ^ 0, where iJ is a static magnetic field, 
X-+(0,0) = xL(O,0) - dM/dH\H=o, where x'.Jq,^) 
is the real part of the longitudinal component of the dy- 
namic spin susceptibility and the magnetisation, M is 
given by2^ 

M = ^E((^kT>-(nki)), (14) 

k 

where nko- is the usual number operator. Writing nto- in 
the spectral representation and performing the sum over 
k one finds that 

^'^=^E'^y_ l^pA^mnpAE), (15) 



where we take cr = ±1, pa{E) = X^k ^so- (k, i?) is the 
density of states of spin a electrons and A^a (k, E) is the 
spectral function for spin a electrons. Even for small 
magnetic fields the Fermi surface may, in general, be 
distorted^^ as E^,^ = E'k - o-i?r(k, Sk^; k, £;k<T)^7e/2 + 
0{H^). However, in the absence of vertex corrections 
such complications cannot arise as the k-dependence 
drops out of the above equation. Therefore we deal with 
the magnetic field by introducing a spin dependent chem- 
ical potential: ^x^ = fi — aHfije/'^ + 0{H^) and noting 
thatriFcr — {l + exp P[{E — ij,cr)]}~^ ■ Thus one finds that. 



hm M = — > cr / — 



~p.{E)Q(E - + a^~p,[E)5{E - ^i„)H + ^^§^Q{E - M.) + 0{H') 



.(16) 



In the limit if — > the first term vanishes due to spin symmetry and the third term vanishes from the assumption 
that we are at an energy scale on which the density of states is constant. Thus, the Knight shift is 

A-..^^, (17) 

which is independent of temperature. 

It follows immediately from Eqs. (llSp and ([T7|) . that JC = 1 for interacting electrons with a contact hyperfine 
coupling when T <C Tq, Tp and neglecting vertex corrections. We stress two points about this result: firstly, our result 
includes the special case of Korringa's result for the free electron gas, /Cftoc = 1; secondly, and more importantly, the 
Korringa ratio is unity for a broad class of systems and not just for the free electron gas. 

Thus we have shown that any deviation of the Korringa ratio from unity, or temperature dependence of 1/TiT or 
Kg must be caused by either vertex corrections, the wavevector dependence of the hyperfine coupling, or the fact the 
there is no energy, fcsTo, on which we may treat the density of states as constant. Note that this result is independent 
of the dimensionality of the system. Further, this is true for any form of the self energy. However, it should be noted 
that the self-energy and the vertex function are not really independent: they both arise from the same underlying 
interactions and can often be related by Ward identities. So this last statement should be taken with appropriate 
caution. 

Let us now briefly discuss the role of vertex corrections in strongly correlated systems. A simple and a widely 
studied approach to nearly (anti)ferromagetic metals is the random phase approximation (RPA). For the Hubbard 
model the RPA givesj^ 




where U (dashed lines) is the effective Coulomb interac- 
tion between electrons on the same lattice site, xo is the 
dynamic susceptibility in the absence of vertex correc- 
tions, and we have suppressed momentum labels in the 
diagrams for clarity. The distinctive 'Stoner-like' form 
of Eq. arises from the sum of ladder diagrams, i.e., 
vertex corrections. For this form of xi^j '-^) ^^'^ ^ three- 



dimensional parabolic band the Korringa ratio is always 
less than unity and decreases monotonically towards zero 
as the Stoner instability is approached. However, for 
a two-dimensional parabolic band the free electron spin 
susceptibility is momentum independent. Consequently, 
c.f. Eqs. (HI), the Korringa ratio remains unity for all U 
even though 1/Ti and the Knight shift both diverge as 
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the Stoner instability is approached.^ 

Korringa type relations have been derived for the im- 
purity spin susceptibility, Ximp{t^), of the Anderson single 
impurity model. It is found tha t^^i^^ 



lim 



It 



(19) 



which is sometimes referred to as the Shiba relation. This 
relation holds even though there can be significant vertex 
corrections for the spin susceptibility. This and the case 
of the two-dimensional RPA with a parabolic band illus- 
trate that the absence of vertex corrections is sufficient 
but not necessary for the Korringa ratio to be unity. 

A similar form of x(q, oj) to that given in (|18p is found 
in dynamical mean field theory (DMFT), but with U re- 
placed by a self-consistently determined four-point ver- 
tex function. In DMFT, as in the RPA, this functional 
form results from a sum over ladder diagrams, with the 
four-point vertex function now forming the legs of the 
ladders. DMFT for the Hubbard model is equivalent 
to an Anderson single impurity model which is solved 
self-consistentlyi^ It follows that the local spin suscep- 
tibility xiocii^) = SqX(q:'^) uiust Satisfy the Shiba re- 
lation ITOl) with Ximp{^) = Xioc{^)- Since, in general, 
Xioci'-^) 7^ x(q = 0,w) the Korringa relation does not 
necessarily hold in DMFT [see Eq. (73) in Ref. 

A slave-boson large N treatment of the periodic An- 
derson model shows that in the low-temperature heavy 
fermion phase the Korringa ratio is close to unity, even 
though 1/TiT and Kg can be several orders of magnitude 
larger than the value predicted in the absence of strong 
electronic correlations ^22 

In the singlet superconducting states the opening of 
the gap causes a suppression of both 1/TiT and Kg- 
For a fully gapped superconductor both 1/TiT and Ks 
show activated behaviours, while if the order parame- 
ter has nodes power laws are seen at low T.^^ The be- 
haviour, particularly of the Knight shift, is more com- 



plicated in a triplet superconductor In the same way, 
a (pseudo)gap destroys the energy scale Tq and hence 
causes non-Korringa behaviour by suppressing 1/TiT 
and Ks- 

Given the above discussion it is interesting to consider 
the normal states of a few superconductors which display 
non-Korringa behaviour. A strong temperature depen- 
dence and a Korringa ratio significantly larger than unity 
is observed in the organic charge transfer saltSji^i^S e.g., 
in K-(ET)2Cu[N(CN)2]Br (Refs. MM)- Heavy fermion 
compounds^ and cobaltates^^ also show enhanced Kor- 
ringa ratios above their Kondo temperatures. Most re- 
cently non-Korringa behaviour has been observed in the 
iron pnictidesi^ In the cuprates, the Korringa ratio of 
YBa2Cu3 6.64 (Refs. US 

has a strong temperature 
dependences and is larger than unity. However, in the 
cuprates the hyperfine coupling has a significant wave 
vector dependencCfS^i^ This complicates deducing the 
relative importance of vertex corrections. We also note 
that vertex corrections are required in the superconduct- 
ing state in order to preserve gauge invariance.'^ 

Moriya's self-consistent renormalised theory^** and the 
phenomenological MMP theoryi^ give a good description 
of many features of the cuprateSfi^i^ii^ organics;^? and 
heavy fermion materials)^. These theories posit a form 
of x(q, oj) which follows from the form of Eq. (|18p . Hence, 
these theories implicitly include vertex corrections. 

In summary, we have shown that for a system with 
a contact hyperfine coupling the Korringa ratio is unity 
in the absence of vertex corrections provided there is an 
energy scale on which the density of states may be treated 
as constant. At sufficiently low temperatures 1/TiT and 
the Knight shift are both independent of temperature 
under the same assumptions. 
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